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Abstract. We present an implementation of the model of minimal universal extra 
dimensions (MUED) in CalcHEP/CompHEP. We include aU level- 1 and level-2 Kaluza-Klein 
(KK) particles outside the Higgs sector. The mass spectrum is automatically calculated at 
one loop in terms of the two input parameters in MUED: the inverse radius R^^ of the 
extra dimension and the cut-off scale of the model A. We implement both the KK number 
conserving and the KK number violating interactions of the KK particles. We also account 
for the proper running of the gauge coupling constants above the clectrowcak scale. The 
implementation has been extensively cross-checked against known analytical results in the 
literature and numerical results from other programs. Our files are publicly available and 
can be used to perform various automated calculations within the MUED model. 
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1. Introduction 

The Standard Model (SM) of particle physics has been successfully verified by experiment 
at low energies. Nevertheless, even if the Higgs boson is discovered, the SM will still be 
considered to be an incomplete theory, as it fails to provide the long-sought missing link 
between Einstein's General Relativity and Quantum Mechanics. The leading candidate for 
a quantum theory of gravity, string theory, typically posits the existence of several new 
ingredients, which are absent in the SM: new spatial dimensions, a symmetry between 
bosons and fermions (supersymmetry), as well as new gauge interactions. All of these new 
ingredients are manifestly present at the Planck scale, but it is not at all clear which of them 
survive down to low energies. Traditionally, supersymmetry and extra gauge interactions 
have attracted the most attention, and their consequences for collider phenomenology have 
been extensively studied [HE]- Within the last 10 years or so, there has been a resurgence of 
interest in models with extra spatial dimensions, whose presence might be revealed in high 
energy collider experiments such as the Tevatron at Fermilab, the Large Hadron Collider 
(LHC) at CERN, or the proposed International Linear Collider (ILC). By now a whole 
plethora of extra-dimensional models have been described and studied to various extent in 
the literature. Roughly speaking, they can all be classified according to the following two 
criteria: 

• How many and which of the SM particles can access the extra dimensions (the bulk). 
The two extremes here are provided by the "large" extra dimension models (also known 
as ADD, after the initials of their original proponents) [3], in which only gravity can 
enter into the bulk, and the Universal Extra Dimensions (UED) models [4], in which all 
SM particles are allowed to propagate in the bulk. 

• What is the metric of the bulk. It can be flat (e.g. in UED), or warped [5]. 

In this paper, we shall concentrate on the simplest case of a single fiat extra dimension, 
which is accessible to the full SM particle content |4] (see Refs. P4TT] for the case of two 
universal extra dimensions). This particular scenario has recently been studied in relation to 
collider phenomenology [T21429] , indirect low-energy constraints [30lfl3] , dark matter [HHSH] 
and cosmology [67H73] . It is therefore of great interest to have an implementation of the 
Minimal UED model (reviewed below in Section [2]) in the most popular general purpose 
computer programs for collider and astroparticle phenomenology. The main goal of this 
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paper is to present one such implementation, suitable for either CalcHEP [71] or CompHEP [75] . 
There are several advantages of choosing CalcHEP and CompHEP for this purpose: 

• CalcHEP and CompHEP can be used for parton- level event generation, preserving the full 
spin correlations in both production and decay. 

• CalcHEP and CompHEP can be easily interfaced [TB] to a general purpose event generator 
such as PYTHIA [77] for the simulation of fragmentation, hadronization and showering. 

• CalcHEP and CompHEP can be easily interfaced with a dark matter program such as 
micrOMEGAs [78] for the calculation of the relic density and detection rates of a generic 
dark matter candidate. 

• The implementation of new models is very straightforward and user- friendly, as we shall 
demonstrate below with the example of Minimal UED. 

The paper is organized as follows. In Section [2] we first review the Minimal UED model 
(MUED), introducing the relevant new particles, couplings and interactions. In Section [3] 
we explain how those were incorporated in CalcHEP and CompHEP. Throughout the paper we 
assume that the readers are already familiar with these programs, so that we only need to 
explain the additional *.mdl model files related to our UED implementatioJll. In Section H] 
we discuss how the implementation can be used to study the collider phenomenology of 
MUED and show some illustrative results. In the Appendices we list some more technical 



results which may be useful to some readers. For example. Appendix A contains the five- 



dimensional UED Lagrangian and Appendix C contains the resulting Feynman rules for the 
level 1 KK particles after compactification. 



2. The Minimal UED Model 

2.1. KK decomposition 

The five-dimensional (5D) UED model [1] is simply the Standard Model placed in an extra 

dimension compactified on an Si/Z^ orbifold, as shown in Fig. [H Let us label the usual 

3 + 1 space-time dimensions with x^, /i = 0, 1, 2, 3, reserving the coordinate y for the extra 

dimension. In order to end up with chiral fermions in 4 dimensions and to project out 

I Our implementation was originally developed for the Second MC4BSM workshop in Princeton, March 
24-27, 2007. Since then, the Minimal UED model has been partially implemented in PYTHIA [79], and more 
fully in CalcHEP, MadGraph, PYTHIA or Sherpa through FeynRules [SOllST]. 
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Figure 1. The S1/Z2 compactification of a single extra dimension on a circle with opposite 
points identified, as indicated by the grey arrows. The blue dots represent the fixed 
(boundary) points and y is the coordinate along the extra dimension. 

unwanted gauge degrees of freedom, one typically imposes an additional symmetry, thus 
creating a manifold with boundaries. For example, in the case of the SijZi orbifold shown 
in Fig. [H one identifies the opposite points on the circle, which creates two fixed points, 
denoted with the blue dots. Any 5-dimensional field can now be assigned a definite parity 
with respect to the orbifold projection : y ^ —y. For example, consider a generic scalar 
field (j){x,y). An even scalar field (j)'^{x,y) is expanded in Kaluza-Klein (KK) modes as 

1 2 °° ny 

(P^ix,y) = ^=(f)^{x) + ^=Y,<Ptix)cos — , (1) 

and obeys Neumann boundary conditions at the two fixed points: 
fd(t)+{x,y)\ fd(t)+{x,y)\ 



dy V 9y 



0. (2) 



y=0 \ y / y=ivR 

Here x is the usual 4-dimensional spacetime coordinate x^, Ris the size of the extra dimension 
and n labels the KK-level. The SM modes correspond to n = 0. In contrast, the KK 
decomposition of an odd scalar field 

= ^= £c(2;)sin^, (3) 

is missing a zero mode (n = 0) and obeys Dirichlet boundary conditions 

0-(x,O) = 0~(x,7ri?) = 0. (4) 

One can similarly assign a definite V5 parity to each component of a gauge field Am{x, y), 
M = 0,1, 2, 3, 5. The usual 3 + 1 components A^, /i = 0, 1, 2, 3, are chosen to be even, which 
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ensures the presence of the SM gauge fields A^^{x) at the n = level, while the extra- 
dimensional component is taken to be odd. The corresponding KK expansions of the 
5- dimensional gauge fields are given by 

A,{x,y) = -^S^A'^{x) + V2f:A;{x)cos{^)^ , (5) 



Mx,y) = ^j:^tAUx)sin{^). (6) 

At the two fixed points y = and y = ttR, the components A^{x, y) (^5(0;, y)) obey Neumann 
(Dirichlet) boundary conditions analogous to eq. ([2]) (eq. (jl])). 

The KK decomposition of a fermion is rather interesting. Since there is no chirality 
in 5 dimensions, the KK modes of the SM fermions come in vectorlike pairs, i.e. there is 
a left-handed and a right-handed KK mode for each SM chiral fermion. For example, the 
5'f/(2)iy-singlet chiral fermions ip%{x) of the SM (which happen to be all right-handed) are 
obtained from the following decomposition 

y) = ^==^l{x) + ^ E rAx) cos ^ , (7) 

rn{x,y) = ^f:rL{x)sm'^, (8) 
V ttR „=i R 

where upon compactification, the two KK fermions ip^x) and ip2ix) given KK level 

n pair up to give a Dirac fermion of mass ^. Similarly, the S'?7(2)vy-doublet SM fermions 
(x) (which happen to be left-handed) arise from 

nix, y) = ^^l(^) + 4^ e n{x) cos ^ , (9) 

*Z(a:,y) = ^E*Msinf , (10) 

where the massive Dirac fermion at each n is now formed from '^^{x) and \l/^(x). 

From eqs. (iTlfTOl) we see that there exist left-handed KK modes ipi{x), which are 
associated with the right-handed SM fermions ipfi{x) and vice versa — there are right- 
handed KK modes \l/^(x), which go along with the left-handed SM fermions (x). This 
often leads to some confusion in the literature when it comes to the labelling of fermion KK 
partners. It should be understood that the chiral index (L or R) of a KK mode fermion 
refers to the chirality of its SM partner. Here we shall also utilize an alternative convention, 
introduced in [82], where the KK fermions are identified by their SU(2)w quantum numbers 
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Table 1. Fcrmion content of the Minimal UED model. S'[/(2)vv'-doublets (S'[/(2)vi/-singlets) 
are denoted with capital (lowercase) letters. KK modes carry a KK index n. and for 
simplicity we omit the index "0" for the SM zero modes. 
SU{2)w representations SM mode KK modes 

. Ul(x) \ ~ ( mix) \ ~ ( mix) 
Quark doublet q^x) = \ ^' Ql{x) ^ \ \ , Ql{x) = \ 

Dl{x) j V Dl{x) j ^ \ Dlix) 

Lepton doublet Ll(x) = Ljix) = \ , L„(x) = 

' El{x) j \ Elix) j " \ Elix) 

Quark Singlet ur{x) Ur{x), u^ix) 

Quark Singlet djiXx) d^x), d^{x) 

Lepton Singlet eii{x) ^r{x), e2(x) 



instead: S'f/(2)vi/-doublets (5'f/(2)i^-singlets) are denoted with capital (lowercase) letters. 
This convention was already employed in eqs. (ITl fTU]) as well. With those conventions, the 
fermion content of the Minimal UED model is listed in Table [H 

Finally, notice that the geometry in Fig. [T]is still invariant under the interchange of the 
two fixed points. The corresponding Z2 symmetry is the celebrated KK parity and will be a 
symmetry of the Lagrangian as long as one continues to treat the two boundary points in a 
symmetric fashion. 

2.2. KK mass spectrum 

At tree level, the mass m„ of any KK mode at the n-th KK level is given by 

^1 = ^ + ^0^ (11) 

where R is the radius of the extra dimension as illustrated in Fig. [H and mo is the mass 
of the corresponding SM particle (zero mode). The resulting mass spectrum for the first 
KK level is shown in Fig. [JK for = 500 GeV, and can be seen to be highly degenerate. 
In fact, several of the lightest n = 1 KK modes have no allowed decays and are absolutely 
stable. 

However, this drastic conclusion is completely reversed, once radiative corrections are 
taken into account [82] • First, the mass spectrum gets renormalized by bulk interactions, 
which are uniquely fixed in terms of the SM gauge and Yukawa couplings, and thus contain no 
new parameters beyond those already appearing in the SM. At the same time, the KK masses 
also receive contributions from terms localized on the boundary points (the two blue dots in 
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Figure 2. The spectrum of the first KK level at (a) tree level and (b) one-loop, for 
R^^ = 500 GeV, AR = 20, = 120 GeV, and assuming vanishing boundary terms at the 
cut-off scale A. (From Ref. [H].) 



Fig. [T]). The coefficients of the boundary terms are in principle new free parameters of the 
theory. The Minimal UED model makes the ansatz that all boundary terms simultaneously 
vanish at some high scale A > R^^. The boundary terms are then regenerated at lower scales 
through RGE running, and lead to additional corrections to the KK mass spectrum [82]. The 
resulting one-loop corrected mass spectrum is shown in Fig. The mass splittings among 
the different n = 1 KK modes are now sufficiently large to allow prompt cascade decays to 
the lightest KK particle (LKP). For the parameter values shown in the figure, the LKP turns 
out to bq§|the KK "photon" 71, although at larger mh the LKP can also be the charged KK 
Higgs boson iff [23] . 

The mass eigenstates of the KK photon 7„ and the KK Z-boson Zn are mixtures of the 
corresponding interaction eigenstates: the KK mode i?„ of the hypercharge gauge boson and 
the KK mode of the neutral SU{2)^r gauge boson. The mixing angle On is obtained by 
diagonalizing the mass matrix in the Wj^) basis 



V 



i + \9lv^ + ^^m^ 



(12) 



where gi ((72) is the hypercharge (weak) gauge coupling, v = 246 GeV is the vev of the SM 
Higgs boson, and 6 represents the total one-loop correction, including both bulk (5) and 

§ Strictly speaking, the true LKP in Fig. ^jp is the KK graviton Qi (not shown). However, due to its 
extremely weak couplings, Qi is irrelevant for collider phenomenology. For its astrophysical implications, 
see [83l. 
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boundary (6) contributions [82]: 

5ml^ = bm\^ + bm\^ . (13) 

Note that for n > 1 the KK mixing angle is in general different from the zero-mode 
(Weinberg) angle 6*0 = dw in the SM. For typical values of and A, 6'„ <^ 6'vi/, and 
the neutral gauge boson KK mass eigenstates become approximately aligned with the 
corresponding interaction eigenstates: 7„ ~ -B„ and Z„ Wj^ for n > 1. This approximation 
will be used in our MUED implementation described below in Section EJ 



2.3. KK interactions 

The bulk interactions of the KK modes are already fixed by the SM. The 5D MUED 
Lagrangian is a straightforward generalization of the SM Lagrangian to 5 dimensions, as 



discussed in [Appendix A[ Upon compactification, integrating over the extra-dimensional 



coordinate y, one recovers the bulk interactions among the various KK modes and their SM 



counterparts (see Appendix C). Since translational invariance holds in the bulk, all these 
bulk interactions conserve both KK number and KK parity. 

However, as already alluded to in the previous subsection, there may also exist 
"boundary" interactions localized on the fixed points in Fig. [H They do not respect 
translational invariance and therefore break KK number by even units. Such interactions 
may already appear at the scale A, being generated by the new physics which is the ultraviolet 
completion of UED. In the Minimal UED version, one makes the assumption that no such 
terms are present at the scale A. Even so, upon renormalization to lower energy scales, 
boundary terms are radiatively generated from bulk interactions. This is illustrated in Fig. |3l 
where we show how an effective coupling between a level 2 KK gauge boson V2 and two SM 
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Table 2. Boundary interactions involving level 2 KK gauge bosons and two SM fermions. 
Here I3 is the fermion isospin and Yt {Yr) is the hypercharge of a left-handed (right-handed) 
SM fermion. In the case of top quarks, one has to include in S(rnf^) the additional corrections 
proportional to the top Yukawa coupling ht: dhtiriT,^ and Sh^mt^, respectively (see [82| for 
details) . 

n ~ 2 KK boson rt = SM fermion Vertex 
U{1)y gauge boson Lepton ^gl7-75T6W 1^ (^)' [^Pl (§5? + f 5I) 

f'^pn m)] 

B2 Quark (up) jg,^M^_L_ in (a)' [^p^ (Xg2 + ^^2 + gg^) 

+ ^PB.{f9! + G9l)] 
Quark (down) ^g,^^^^^\n (a)' [^p^ (r_g2 + + gg^) 

+^^i^(|g? + 6gi)] 

5f/(2)w/ gauge boson Lepton 1/3327^75 leW In (^) ^'l [fff? - f 5I] 





Z2 


Quark 


*/3327-^TaWln(^) 


'PL[yi-f9'2+G9l] 


5[/(2)h. 


gauge boson 
W2 


Lepton 
Quark 


^^7-^TeWln(^)' 


Pl [hi - f 5i] 
PL[yi-f9l + G9l] 


5[/(3)e 


gauge boson 


Quark (up) 




)'[PL(ig? + f.9i-¥5l) 



G2 Quark (down) ^53^y'-l=^ In (a)' [p^ (ig2 + 27^2 n^2) 



fermions is generated at one loop from a diagram with level 1 KK particles running in the 
loop. This effective coupling 

a/2 \ mi m2 J 

can be expressed in terms of the boundary contributions 5m„ (see eq. ( IT5]) ) to the one- loop 
mass corrections [82] • The explicit form of this effective coupling is summarized in Table [2] 
for each different type of level 2 KK gauge boson and for the various possible SM fermion 
pairs. 

3. Model files 

Having reviewed the MUED model, we are now in a position to describe its implementation 
in CalcHEP and CompHEP. Each one of these programs gives its users an opportunity to 
incorporate new physics in the already existing framework of the SM, MSSM, etc. To 
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Table 3. KK gauge bosons. 



Name 


A 


A+ 


2*spin 


mass 


width 


eolor 




KG 


KG 


2 


MKG 


wKG 


8 




Bl 


Bl 


2 


MBl 





1 




Zl 


Zl 


2 


MZl 


wZl 


1 




- W+ 


^ M^- 


2 


MWl 


wWl 


1 


Gl 


~ G2 


~ G2 


2 


MKG2 


wKG2 


8 




B2 


B2 


2 


MB2 


wB2 


1 


zl 


Z2 


Z2 


2 


MZ2 


wZ2 


1 


wl 


- W2 


- w2 


2 


MW2 


wW2 


1 



this end, one must simply supply an updated version of the four model files defining a 
given physics scenario in CalcHEP and CompHEP: prtclsN.mdl, varsN.mdl, funcN.mdl and 
IgrngN.mdl, where N stands for the numerical label of the physics scenario in the model 
menu of CalcHEP and CompHEP. We shall now discuss each one of those files, which are 
available from [http : //home . f nal . gov/ ^kckong/mued/ 

3.1. Particles 

New particles are defined in the prtclsN.mdl model file. We incorporate the n = 1 and 
n = 2 KK modes of the gauge bosons (see Table [3]), leptons (see Table H]) and quarks (see 
Table [5]). In Tables [3l|5] the KK number is represented by a superscript n = 1 or n = 2, while 
the subscript is either the Lorentz index (/i) of the vector particles in Table [3] or the chirality 
index of the fermion particles in Tables H] and [5l We remind the reader that all KK fermions 
are vectorlike and the chirality index refers to the chirality of their SM counterparts. The 
corresponding masses and widths of the KK fermions in Tables H] and |5] carry "D" or "S" to 
indicate their nature, SU {2)\Y-donhlet or S'?7(2)vi/-singlet, respectively. The new particles 
listed in Tables [3ll5] are in addition to the usual SM particles which are not shown here. 

3.2. Variables 

The input parameters for any given physics scenario are defined in the varsN.mdl model file. 
In principle, MUED has only two additional input parameters beyond the SM: the radius R 
of the extra dimension and the cut-off scale A. For convenience, we use the inverse radius 
and the number of KK levels AR which can fit below the scale A. R~^ has dimensions 
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Table 4. KK leptons. 



Name 


A 


A+ 


2*spin 


mass 


width 


color 


ei 


- eL 


- EL 


1 


DMe 


wDel 


1 


Mi 


~ mi 


- ML 


1 


DMm 


wDe2 


1 


rh 


- tL 


- TL 


1 


DMt 


wDe3 


1 


ek 


- ei? 


- Li? 


1 


SMe 


wSel 


1 


Mi? 


^ mi? 


- MR 


1 


SMm 


wSe2 


1 




- m 


- Ti? 


1 


SMt 


wSc3 


1 




nl 


- 7V1 


1 


DMcn 


wDnl 


1 




- n2 


- N2 


1 


DMmn 


wDn2 


1 




n3 


- TVS 


1 


DMtn 


wDn3 


1 


el 


~k 


- IE 


1 


DMe2 


wDel2 


1 


t^l 


'--^ Im 


- /Ai 


1 


DMm2 


wDe22 


1 


rl 






1 


DMt2 


wDc32 


1 


el 


~ re 


- rL 


1 


SMe2 


wScl2 


1 




^ rm 


- rM 


1 


SMm2 


wSe22 


1 


rl 


- ri 


- rT 


1 


SMt2 


wSe32 


1 


-1 


^ en 


- eTV 


1 


DMen2 


wDnl2 


1 


-I 


mn 


- miV 


1 


DMmn2 


wDn22 


1 




^ tn 


- tN 


1 


DMtn2 


wDn32 


1 



of GeV, while AR is dimensionless. Our additions to the varsN.mdl model file are listed in 
Table ini As seen from the table, we also include several other variables of interest. RG is used 
to turn on and off the running of coupling constants, while scaleN is the renormalization 
scale fi at which the couplings are evaluated. The remaining parameters in Table |6] are 
some useful numerical constants related to the RGE running of the gauge couplings (see 
Section [3751) . 

3. 3. Constraints 

The funcN.mdl model file is reserved for variables which are not numerical inputs, but 
are instead computed in terms of the parameters already defined in the varsN.mdl model 
file. In our case, we use funcN.mdl to supply the masses and two-body decay widths of 
the KK particles introduced in Section 13.11 Therefore they are automatically computed 
by CalcHEP/CompHEP at the beginning of each numerical session. The masses for all KK 
particles are evaluated based on the 1-loop formulas of Ref. [H2] and we have also made 
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Table 5. KK quarks. 



Name 


A 


A+ 


2*spin 


mass 


width 


color 




Du 


DU 


1 


DMu 


wDu 


3 




Dd 


DD 


1 


DMd 


wDd 


3 




Dc 


DC 


1 


DMc 


wDc 


3 


s\ 


Ds 


DS 


1 


DMs 


wDs 


3 


t\ 


Dt 


DT 


1 


DMtop 


wDt 


3 




Db 


DB 


1 


DMb 


wDb 


3 




Su 


SU 


1 


SMu 


wSu 


3 




Sd 


SD 


1 


SMd 


wSd 


3 




Sc 


SC 


1 


SMc 


wSc 


3 


4 


Ss 


SS 


1 


SMs 


wSs 


3 




St 


ST 


1 


SMtop 


wSt 


3 




Sb 


SB 


1 


SMb 


wSb 


3 




- Du 


- DU 


1 


DMu2 


wDu2 


3 




- Dd 


- DD 


1 


DMd2 


wDd2 


3 




- L>c 


- DC 


1 


DMc2 


wDc2 


3 


s2 


- Ds 


- DS* 


1 


DMs2 


wDs2 


3 




- Dt 


- DT 


1 


DMtop2 


wDt2 


3 




- £16 


- DB 


1 


DMb2 


wDb2 


3 


Ur 


~ Su 


- s-c/ 


1 


SMu2 


wSu2 


3 


dl 


- 5d 




1 


SMd2 


wSd2 


3 




- S-c 


- s-c 


1 


SMc2 


wSc2 


3 




- 5s 




1 


SMs2 


wSs2 


3 




- Si 


- ST 


1 


SMtop2 


wSt2 


3 


^1 


- Sh 


- 55 


1 


SMb2 


wSb2 


3 



numerical cross-checks with the results from the private code used in Ref. [82]. Our formulas 
for the widths have been derived analytically and cross-checked with CalcHEP/CompHEP (see 
Section H]). A partial list of 2 body decay widths can be found in p ^ [T5 | [20] and our formulas 
agree with their expressions. In the older versions of CalcHEP/CompHEP, defining the widths 
as constraints was very convenient in our implementation, since one did not have to launch 
a separate numerical session for their calculation, and then enter their numerical values as 
input parameters. However, the more recent versions of CalcHEP and CompHEP allow for the 
automatic calculation of the particle widths on the fly, using the interactions defined in the 
IgrngN.mdl model file. Our implementation thus allows for backward compatibility with 
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Table 6. Parameters added to the varsN.mdl model file. 



Parameters 


Default values 


Symbols 


Comments 


Kmv 


500 


D-l 


Inverse radius of the extra dimension 


LR 


20 


AR 


The number of KK levels below A 


RG 


1 




1 turn on the running of the coupling constants 
U turn on the running of the coupimg constants 








Renormalization scale, ^ = 


scaleN 


2 


n 


n=2 can be used for KK level 1 pan- production 

or level 2 single production 
11=4 can be used for KK level 2 pair production 


cbl 


6.8333 




^, The coefficient of the SM /5-function for U[1)y 


cb2 


-3.16667 




— The coefficient of the SM /3-function for 5f7(2)vv 


cb3 


-7 




-7, The coefficient of the SM ;3-function for SU[3)c 


cblt 


6.8333 




The coefficient of the KK /3-fuiiction for U{1)y 


cb2t 


-2.83333 




— The coefficient of the KK p-function for SU[2)w 


cb3t 


-6.5 


h 


The coefficient of the KK /3-function for SU{3)c 


clMZ 


98.4151 






c2MZ 


29.5846 






c3MZ 


8.53244 




a^\pi^Mz) 



older versions of CalcHEP/CompHEP. 
3.4- Interactions 

The new interactions of the KK particles of Section 13.11 are added to the IgrngN . mdl model 
file. We include the usual bulk interactions, as well as the KK number violating boundary 
interactions listed in Table [2] [82] • Since the Weinberg angle On for any n > 1 is small [82] , 
we ignore the mixing among the neutral KK gauge bosons. Thus the KK-photon 7^ is 
identical to the hypercharge gauge boson i?„ and the KK Z-boson is identical to the 
neutral SU{2)w gauge boson W^- We also ignore the mixing between S't/(2)vK-doublet and 
S'f/(2)vi/-singlet KK fermions. 

Our IgrngN. mdl model file includes all interactions of level-1 and level-2 KK particles 
except for the KK Higgs bosons. The phenomenology of the KK Higgs bosons is very model 
dependent, depending on the value of the SM Higgs mass rrth and the bulk Higgs mass term 
(see |82] for details). Therefore we omit any interactions involving KK Higgs bosons 
II The collider phenomenology of the KK Higgs bosons has been discussed in [23 1 124 1 129]. 
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The UED Lagrangian can be easily derived as shown in [Appendix C[ Here we only 



point out how to deal with 4-point interactions involving KK gluons, since this case requires 

special treatment when implemented in CalcHEP/CompHEP. 

The Lagrangian for the quartic interactions with KK gluons is the following 
1 

n _ 2 rabc radeQO^bQO,CQOjdfi^O,eiy ^3 j^abc rade^O,d^O,e^l,bfi^l,ch' 

- f {r'^G'^GY + G'fGf))' - ^ ■ ^-glr'^^r'^^G'^GYG'^'^G'^^'^ . (14) 

The color structure of these 4-point interactions cannot be directly written down in 
CalcHEP/CompHEP format. Hence, to implement this vertex in CalcHEP/CompHEP, we use 
the following trick. We introduce three auxiliary tensor fields t^^,, s^^ and in the same 
way as the original CalcHEP/CompHEP approach for SM gluons. Then one can rewrite the 
Lagrangian as 

C = - -eX^"" + i^t':,J''^^G^^''G^^'' + i-^tU'^^^G^^^'G^''' 
- -sLs""^ + i—s'^^J^^^G^^^'G^^" 



-u,„,u 



2 '72 



- \ \t% - ^g^^r'^'G'^G': - tg.^r'^G'^Gl'^ (15) 

yabc jade f^Ob^Oc _|_ (^^^ (G^'^^G^^'^ -\- G^'^^ G^^^ 



1 / 1 \2 1 

J- / „ .J- „„h^^^h^^^,\ J- 



4- 

2 ['lu - ^gs^r'^Gl'Gl^j - ^glr'^r'^Gl'Gl^G'^^G'^ 

- 1 - m^r'" {Gl'G'; + Gl'Gl^)^ 

-\9l{r''{Gl'Gl^ + Gl'Gl^))\ 

It is easy to show that the functional integration over the three auxiliary tensor fields 
reproduces the 4-gluon interactions (fT4l) . 

3. 5. Running of the coupling constants 

Due to the additional contributions from the KK modes to the beta functions, the gauge 
couplings run faster in theories with extra dimensions. The RGE for ctj = |^ is given by 



daj ^ _ hj - hj _ hjXs ' ^ ^ 
dt 2-71 2iT \/io. 



(16) 
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where S is the number of extra dimensions, /iq is some reference energy scale, Xs = ^f^^y^, 

/41 19 \ 

{buh,h) = [-,--,-7) (17) 
are the SM beta function coefficients, while 

correspond to the contributions of the Kaluza-Klein states at each massive KK excitation 
level [85l[86]. The solution to ( IT6|) becomes 

-1 \ 1 /i biXs 



\ <5 



-1 . (19) 

The effect of the RGE running f|T9|) can be accounted for by setting the RG parameter in 
Table [6] to 1 and choosing the appropriate renormalization scale via scaleN. 

4. Discussion 

4-1. Code validation 

In general, the availability of CalcHEP/CompHEP model files opens the door to a number 
of applications related to collider phenomenology and dark matter searches. Each such 
individual study contributes to the validation of the code. Further consistency checks are 
provided by comparing to existing analytical and/or numerical results in the literature. 

• For starters, we have compared the KK mass spectrum calculated with our 
implementation to the results shown in Fig. which were obtained independently in 
Ref. [82]. Using identical inputs, and neglecting the running of the gauge couplings (as 
was done in [HS]), we found perfect agreement. 



The interaction vertices of Appendix C can be independently derived with the 



automated tool LanHEP [87]. We checked some of the more technically challenging 
cases (especially the self-interactions of gauge bosons) and also found agreement. 

To minimize the possibility of typing mistakes, we computed analytically the cross- 
sections for a selected number of simple scattering processes, and compared to the 
analytical expressions derived by CalcHEP/CompHEP. 

We have similarly checked that the KK particle widths calculated from our analytical 
expressions agree with those computed with CalcHEP/CompHEP by means of our MUED 
implementation. 



16 




jQ-4 I .. I .... I .... I .... I .... I .... I 

400 600 BOO 1000 1200 1400 
R"' (GeV) 




J 0-5 I ^ \ ^ ^ ^ \ ^ ^ ^ \ ^ ^ ^ \ ^ ^ ^ ^ ^ ^ ^ L 

400 600 800 1000 1200 1400 

R"' (GeV) 



Figure 4. Strong production of ji = 1 KK particles at the LHC for = 7 TcV: (a) KK- 
quark pair production; (b) KK-quark/KK-gluon associated production and KK-gluon pair 
production. The cross sections have been summed over all quark flavors and also include 
charge-conjugated contributions such as QiQi, QiQi, giQi, etc. We use CTEQ6L parton 
distributions |91| and choose the scale of the strong coupling constant as to be equal to the 
parton level center of mass energy. 

• Our analytic formulas for decay widths agree with the expressions given in [T^[T^ [2U]. 

• Our implementation was used for the analytic calculation of all (co) annihilation cross- 
sections of level 1 KK particles (SH] and the results were in complete agreement 
with [111158]. 

• Our model files have already been used for various collider studies p!7 t[T9] - [22| l26 y 59 [ [88]- 
[90] . One example is shown in Fig. IH which shows the strong production cross-section 
of level 1 KK particles at the imminent LHC energy of 7 TeV. 

• We have compared results for various production cross-sections in MUED to those in 
published papers [T21[I2] and find agreement. 

• Our model files were also cross-checked against the known analytical expressions for 
various invariant mass distributions [TSIIH21IH3] • 

• Our model files have also been tested by other groups, for example in creating 
Pythia_UED [79 | l9 ^ l95] , which implemented the matrix elements for certain processes in 
PYTHIA [77]. Another extensive comparison to an independent MUED implementation 
via FeynRules was done in [8T] . 
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4-2. Outlook 

Moving forward, it is important to be mindful of the limitations of our implementation. 
First of all, it is still Minimal UED, and the spectrum is quite constrained, given in terms of 
only 2 parameters: and A. If a signal consistent with UED is discovered at the LHC or 
the Tevatron, one would like to start testing the data with a more general UED framework, 
which allows for the presence of arbitrary boundary terms at the scale A. Work along these 
lines has already started and a beta version of the corresponding UED model files is available 
from the authors upon request. 
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Appendix A. UED Lagrangian in 5 dimensions 
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The Lagrangian for the 5-dimensional UED model is written as 



GF 



ttR 



^Leptons 
^Quarks 



ttR 



Yukawa 



ttR 



ttR 



(A.l) 
(A.2) 



dy{iLix,y)T^'DML{x,y) + ie{x,y)T^'DMeix,y)} , (A.3) 

dy {iQ{x, y)T^'DMQix, y) + iu{x, y)T''DMu{x, y) (A.4) 

+id{x,y)T^^ DMd{x,y)^ , 

dy [\uQ{x, y)u{x, y)iT'^H*{x, y) + \dQ{x, y)d{x, y)H{x, y) 
-XeL{x,y)e{x,y)H{x,y)} , (A.5) 

dy [{DmH{x, y)^ {d^'H{x, y)) + fi^H^{x, y)H{x, y) 

-\{H\x,y)H{x,y)f'^ , (A.6) 

in terms of 5-dimensional fields decomposed as discussed in Section 12.11 



C 



Higgs 



H{x,y) 

B^,{x,y) 

B^{x,y) 

W^{x,y) 

W5{x,y) 

Gf^ix^y) 

G5{x,y) 

Q{x,y) 



y=lH{x) + V2Y.Hn{x)cos{^] 



n=l 

oo 



^lBlix) + V2j:B;{x)cos{'^] 



n=l 



^\w'^{x) + V2Y.W-{x)cos{^^ 



1 



GI{x) + V2Y.G;{x)cos< 



ttR 



EGl 



,Xj sm 



n=l 



ny, 
R' 



n=l 



'txR 



n=l 



,ny. 



(A.7) 



ny. 



PLQl{x)cos{^) + PnQl{x)sin{^ 
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u{x,y) 
d{x,y) 
L{x,y) 
e{x,y) 



1 



^ttR 
1 



n=l L ^ 



1 



1 



k n=l 
' oo r 



'TT 



R 



ny 

r' 

^ny. 
' R' 

eH,(x) + y2^ P^e^(a;)cos(^) 

n=l L ^ 



n=l 

oo 



PRdnix) cos( 
Pz.i^2(^) cos 



-PrL^ (x) sin 



P . 
ny 
^ P 



P,e2(x)sin(|^) 



Here H{x,y) is the 5D Higgs scalar field and {B^{x,y), B^{x,y)), {W^{x,y),W^{x,y)) and 
y), G^{x, y)) axe the 5D gauge fields Bm, Wm and Gm for f/ (l)y, Sf/ (2)vi/ and Sf/ (3)c, 
respectively. The 5D index M runs over M = /i, 5, where = 0, 1,2,3. The SU{2)w and 
S'[/(3)c gauge fields are 



G 



2 



A/" 



where r", a = 1,2, 3, are the usual Pauli matrices and A"^, A = 1, 2, 8, are the usual Gell- 
Mann matrices. The 5D field strength tensors for f/(l)y, SU{2)]^ and SU{3)c are defined 
as follows 



B 



MN 



QmBm — SmB 



w 



MN 



(A.8) 



where e"'^'^ and f'^^'" are the structure constants for SU{2)w and SU{3)c, respectively. The 
parameter ^ in f lA.2p is the gauge fixing parameter in the generalized Pg gauge. 

The 5-dimensional (4-dimensional) gauge couplings are denoted by g^^^ (gi), where 
i = 1,2,3 stands for f/(l)y, SU{2)w and SU{3)c, correspondingly. The two types of 
couplings are related by 



/ttR 



(A.9) 



Finally, Q{x, y) and L(x, y) are the S'?7(2)vy-doublet fermions from Tabled], while u{x, y), 
d{x, y) and e(x, y) are the corresponding S'f/(2)vi/-singlet fermions from Tabled] Pl,r = ^^^^ 
are the 4D chiral projectors in terms of the usual 75 matrix. The gamma matrices in 5D 

^^ = (7^^7'), (A.IO) 
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satisfy the Dirac-Clifford algebra 



{r^,r^} = 2(7 



MN 



where the 5D metric 

/ 

Qmn = 







V 



9^^p 
-1 

and g^'^ = (H ) is the usual 4D metric. 

The covariant derivatives act on 5D fields as follows 



DAiQ{x,y) 
Dmu{x, y) 
DMd{x, y) 
DML{x,y) 
DMe{x,y) 



d 



M 



,(5)l 



Yo 



ig^'Gn + i^gf^BM ] d{x, y) 



.Yi 



dM + + i^-^g?BM ) L{x, y) 

Qm + i^gf^BA.j^ e{x,y) , 



Q{x,y) 



(A.ll) 



(A.12) 



(A.13) 



where the fermion hypercharges are Yg = |, = |, 1^ = — |, Fx, = — 1 and = — 2. 

It is now a rather straightforward but tedious exercise to substitute the expansions (1A.7P 
into the 5D Lagrangians (lA.ltlA.6]) and perform the integration over y with the help of the 



orthonormality relations listed in Appendix B The resulting Feynman rules in terms of 



4-dimensional fields are listed in Appendix C 
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Appendix B. Ortho normality Relations 



The following orthonormality relations can be used in the process of compactifying the 5- 



dimensional Lagrangian listed in Appendix A 

dy cos(— )cos( — 
K K 

^ my . ,ny 







ttR 



rR 



ttR 



ttR 



ttR 



rR 



ttR 



ttR 



dy sin 
dy cos 
dy cos 
dy sin 
dy sin 
dy sin 
dy cos 
dy sin 
dy sin 
dy sin 



rR 



dy sin 



, sm 

my ny ly 
— )cos(— )cos(-) 

my. .ny ly 

1 cosi — I cosi — I COSI — 

r' ^R' ^R' ^R' 

my. . ,ny . ly . ky 

-^) ^^^^W 

my . ny ly 
— )sin(— )cos(-) 

my. . .ny ly ky 
1 sml — I COSI — I cosi — I 

r' ^r' ^r' ^r' 

my . ny 

-R^^'^'^R^ 

my . ny . ly 
— )sin(— )sin(-) 

my ny ly 
_)cos(— )cos(-) 

my ny ly ky 
-^) ) cos(— ) cos(— ) 

my. . ,ny . ly ky 
1 sml — I sml — I COSI — I 

r' ^r' ^r' ^r' 



ttR 

2 "m,n ! 

ttR 

2 "m,n ) 

—A' 

^ mnl > 

o mnlk ' 



ttR 



mnlk ' 



ttR 4 

^ mnl 5 
g mnlk 1 

0, 
0, 
0, 
0, 
0, 



(B.l) 



where the A symbols are defined as 



Al 



mnl 



A2 

mnlk 



A^ 

mnlk 



^mnl 

A^ 

mnlk 



s, 



l,m+n ~l~ ^n,l+m 



+ 6, 



m,l+n ) 



6, 



n + Si 



,m+n+k 



5. 



m,n+k+l 



+ 5n 



,k+l+m 



k+m,l+n + h+i n,l+m 1 

~ Sk,l+m+n ~ Sl,m+n+k ~ Sm,n+k+l ~ Sn,k+l+m 
l+n + 4+ n,l+m 1 
^l,m+n ~l~ ^n,l+m ~l~ ^m,l+n t 
- 6k,l+m+n - Sl^m+n+k + Sm,n+k+l + S n,k+l+m 
~ ^k+l,m+n + Sk+m,l+n + ^k+n,l+m ■ 



(B.2) 
(B.3) 

(B.4) 

(B.5) 
(B.6) 
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Appendix C. Feynman Rules 



Here we list some of the KK-number conserving vertices with subscripts 'n' standing for 
the KK-level, which are obtained after compactifying the 5-dimensional Lagrangian of 
Appendix A| with the help of the ortho normality relations of Appendix B For KK-number 
violating (but still KK-parity conserving) vertices, refer to Fig. |3] and Table [2l 



QC 




ba 



MM' 




-,Sb 



G 




G 



P 



^ ^ = ^93f [{p - q)xgt.u + {q- r)^,gxu + (r - p)^^am] 



Ac 



G 



ub 





G 



Ac 





+ r''f'"''{9x,9.p-9xu9,p) 



GP^ 
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^ 2 V2 ' ' 



ff 



fi 



D 




-ie[{ki - k2)gf,u + (h - h)g^x + (^3 - kl)gx^,] 



W+ 




-ig2[{h - k2)g^,u + {h - k^)g^,x + {kj. - ki)gx^] 




-ig2 cos6'iy[(A;i - k2)g^,u + {k2 - k^)g^x + (/^s - k^ 
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